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ABSTRACT 

Two types of correlations between the radio and X-ray luminosities (Lr and Lx) have 
been found in black hole X-ray binaries. For some sources, they follow the ‘original’ 
type of correlation which is described by a single power-law. Later it was found that 
some other sources follow a different correlation consisting of three power-law branches, 
with each branch having different power-law indexes. In this work, we explain these two 
types of correlation under the coupled accretion-jet model. We attribute the difference 
between these two types of sources to the different value of viscosity parameter a. 
One possible reason for different a is the different configuration of magnetic field in 
the accretion material coming from the companion stars. For the ‘single power-law’ 
sources, their a is high; so their accretion is always in the mode of advection-dominated 
accretion flow (ADAF) for the whole range of X-ray luminosity. For those ‘hybrid 
power-law’ sources, the value of a is small so their accretion mode changes from an 
ADAF to a luminous hot accretion flow, and eventually to two-phase accretion as the 
accretion rate increases. Because the dependence of radiative efficiency on the mass 
accretion rate is different for these three accretion modes, different power-law indexes 
in the Lr - Lx correlation are expected. Constraints on the ratio of the mass loss rate 
into the jet and the mass accretion rate in the accretion flow are obtained, which can 
be tested in future by radiative magnetohydrodynamic numerical simulations of jet 
formation. 

Key words: accretion, accretion discs - black hole physics - ISM: jets and outflows 
- X-rays: binaries 


1 INTRODUCTION 

A strong correlation between the radio luminosity (e.g., at 
8.6 GHz; hereafter Lr) and X-ray luminosity (e.g., at 3-9 
keV; hereafter Lx) has been found in black hole (BH) X-ray 
binaries (BHBs) in their hard state (e.g.. Corbel et al. 2000, 
2003, 2013; Gallo et al. 2003; but see Xue & Cui 2007). 
This correlation can be well described by a single power- 
law form, Lr oc lJ, with the index p ^ 0.6 (Corbel et al. 
2013). We call this the ‘original’ radio/X-ray correlation. 
The correlation was extended to including low-luminosity 
active galactic nuclei (AGNs) by considering the effect of 
the BH mass (e.g., Merloni et al. 2003; Falcke et al. 2004; 
Hording et al. 2006; Wang et al. 2006; Li et al. 2008; Giiltekin 
et al. 2009; Younes et al. 2012). 

However, it was soon found that some sources do not 
follow this single power-law correlation. Rather, they follow 
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a ‘hybrid’ correlation, i.e. the correlation index varies in 
different regimes of luminosity. Among them H1743-322 is a 
prototype, thanks to the frequent outbursts it undergoes and 
the sufficient large coverage in its Lx during the outbursts. 
The radio/X-ray correlation index of this source shows a 
transition from p s; 1.4 at the bright Lx regime, through 
p ~ 0.2 (almost a flat correlation) at moderate Lx regime, 
eventually to the ‘original’ p » 0.6 at the weak Lx regime 
(Coriat et al. 2011). Such kind of hybrid correlation was 
also discovered, but with fewer data points and/or narrower 
coverage in Lx, in MAXI J1659-152 (Jonker et al. 2012) and 
XTE J1752-223 (Ratti et al. 2012; Brocksopp et al. 2013). 
Statistical analysis of 18 BHBs by Gallo et al. (2012) also 
shows that, in addition to the original p x 0.6 correlation, 
another correlation with p x 0.98 exists at the bright Lx 
regime. All these observational results are summarized in 
Fig. 1. 

Theoretically the ‘original’ p x 0.6 correlation has 
been explained in the framework of the coupled accretion- 
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jet model (Heinz & Sunyaev 2003; Merloni et al. 2003; 
Heinz 2004; Yuan & Cui 2005). In this model, the radio 
radiation comes from the synchrotron emission from the jet, 
while X-ray from the Comptonization emission from the 
hot accretion flow. In this case, the correlation index p is 
mainly determined by three factors, namely the dependence 
of Lr on the mass-loss rate in the jet Mjet (i.e., the radio 
radiative efficiency of jet), the dependence of Lx on the mass 
accretion rate Mm (i.e., the X-ray radiative efficiency of the 
accretion flow), and finally the fraction of accreting material 
that enters into the jet (see Eq. 1 below for definition) 
as a function of accretion rate. 

In addition to explaining the observed ‘original’ cor¬ 
relation, two predictions are made by Yuan & Cui (2005) 
based on the accretion-jet model by extrapolating rjja to 
lower accretion rates. One is that when the Lx Lx is below a 
critical value, Lx.crit ~ 10“^ - 10“'’ Ledd (Lsdd is the Eddington 
luminosity), the X-ray radiation from the jet will exceed that 
from the accretion flow. Physically, this is because the jet 
also emits X-ray radiation and this emission is less sensitive 
to the accretion rate compared with the radiation from the 
accretion flow; thus below a very low accretion rate, the 
radiation from the jet will catch up with that from the 
accretion flow and even become a dominant contributor. 
This explains why the X-ray emission of several very low 
luminosity sources such as M87 is dominated by the jet 
(Wilson & Yang 2002; Yuan, Yu & Ho 2009). The second 
prediction is that the radio/X-ray correlation should become 
steeper, i.e. p ~ 1.23 (Yuan & Cui 2005, see also Heinz 2004). 
Both predictions have been confirmed by later works, both 
observational and theoretical (Pellegrini et al. 2007; Wu et 
al. 2007; Pszota et al. 2008; Wrobel et al. 2008; Yuan et 
al. 2009; de Gasperin et al. 2011; Younes et al. 2012). For 
example, it is found that the data from all available low- 
luminosity AGNs satisfying Lx < Lx.crit follows Lr oc in 
perfect agreement with the Yuan & CUi (2005) prediction. 
In the case of BH X-ray binaries, however, the answer is 
less clear due to two reasons. One is that the data points of 
BHBs in quiescent states are still very limited. Another is 
that the data quality of radio observations at such low X- 
ray luminosities is still very poor (Yuan & Narayan 2014). 
For example, the radio detection of XTE J1118+480 at its 
quiescent state is marginal, only at 3cr level (Gallo et al. 
2014). 

In this work, we focus on relatively luminous obser¬ 
vations. The question we want to address is whether we 
can also understand the ‘hybrid’ correlation based on the 
accretion-jet model. Before we begin our study, we would 
like to mention several models proposed recently. Meyer- 
Hofmeister & Meyer (2014, see also Gao, Wu & Dong 2014; 
Huang, Wu & Wang 2014 and Qiao & Liu 2015) propose 
that when accretion rate is high, a weak, cool disc will 
be formed in the innermost region of the hot accretion 
flow. This cool disc supplies additional seed photons for 
Gompton scattering, which will significantly enhance the 
produced X-ray flux. They propose it can explain the p x lA 
radio/X-ray correlation, i.e. the bright part of the ‘hybrid’ 
correlation. However, they do not explain why there is also a 
flat p ~ 0.2 branch and why different sources follow different 
correlations. These questions will be addressed In this work. 

This paper is organized as follows. In Section 2 we 
briefly review the accretion-jet model. We then calculate the 


radiation from both the hot accretion flow and jet to explain 
the correlations in Section 3. The last section (Section 4) 
devotes to a summary and discussion. 


2 THE ACCRETION-JET MODEL 

There are three components in this model, i.e. an outer 
truncated Shakura-Sunyaev disc (Shakura & Sunyaev 1973, 
SSD hereafter), an inner hot accretion flow, and a jet. The 
details of this model are described in Yuan, Cui & Narayan 
(2005, hereafter YCN05). 

2.1 The jet model 

Recently many magnetohydrodynamical (MHD) numerical 
simulations have been performed to study the jet formation 
(see review in Yuan & Narayan 2014). Most of these works 
focus on the ‘BZ-jet’, which is formed by extracting the 
spin energy of the BH [e.g., Blandford & Znajek 1977 (BZ); 
Tchekhovskoy et al. 2011]. It is a relativistic, Poynting 
flux-dominated jet. In addition to the BZ jet, numerical 
simulations have revealed another type of jet called ‘disc- 
jet’, which is formed by extracting the rotation energy of the 
underlying disc. It is a quasi-relativistic, matter-dominated 
jet. The ‘disk-jet’ can even exist around a non-rotating BH 
since it is powered by the rotation of the accretion flow (e.g., 
Lynden-Bell 2003; Kato et al. 2004; Hawley & Krolik 2006; 
Ohsuga & Mineshige 2011; Yuan et al. 2015). 

The jet scenario adopted here seems to be more close 
to the ‘disk-jet’ model, because the composition of the 
jet is assumed to be dominated by normal plasma, i.e. 
electrons and ions, coming from the underlying hot accretion 
flow. The existence of protons in the jet is supported by 
several observations of the AGN jets (see Sikora 2011 for 
a review): (1) the detection of circular polarization and 
Faraday rotation of the radio core. This is because the 
electron-positron plasma cannot generate any polarization 
(Park & Blackman 2010); (2) the low-energy cutoff in the 
radio spectra (and also the electron energy distribution) 
of hotspots in radio-lobes, which is likely a consequence 
of dissipation of bulk kinetic energy in an electron-proton 
jet (Godfrey et al. 2009). We mainly follow Spada et al. 
(2001) to calculate the radiation from the jet*. The half 
opening angle of the jet is assumed to be 0 = 0.1, and 
the bulk Lorentz factor is fixed to Pjc = 1.2. Within the 
jet, internal shocks occur due to the collision of shells with 
different velocities. These shocks accelerate a fraction of the 
electrons = 0.01) into a power-law energy distribution. 
Shock acceleration theory predicts that the power-law index 
Pjc, of these non-thermal electrons to be 2 < pje, < 3, 
and we set pjc = 2.14 throughout this work. The steady- 
state energy distribution of the accelerated electrons is 
then self-consistently calculated, taking into account the 
effect of radiative cooling. The energy density of accelerated 
electrons and amplified magnetic field is determined by two 
parameters, 6c and eg, which describe the fraction of the 
shock energy that goes into electrons and magnetic fields, 

' Recently Kumar & Crumley (2015) calculate of the radiation 
from the ‘BZ-jet’. 
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respectively. Theoretical study of relativistic collisionless 
electron-proton plasma puts a constrain on and eg, i.e. 
fc ~ (Medvedev 2006). We thus take the values of 
and eg to be 0.1 and 0.02, respectively. These values are 
also within the typical range obtained in the study of GRB 
afterglows (Medvedev 2006). We note that different values 
of Cf and eg will only change the normalization, but not 

the slope, of the T/jet-Min(5f?s) relationship (cf. equation 

3 below). We then calculate the synchrotron emission from 
these accelerated electrons. 

Like any other jet models in literature, our jet model 
is phenomenological. Besides, we assume that all the jet 
parameters remain unchanged for different sources and 
accretion rates, which is certainly a strong assumption. This 
is partly because we still do not have good constraint on 
them. To conclude, the only free parameter in our jet model 
is its mass-loss rate Mjet. We note that the Lr is found to 
be a power-law function of (see also Heinz & Sunyaev 
2003). 


2.2 Hot accretion flow model 

Depending on the accretion rate, two types of hot accretion 
flow exist. They are the advection-dominated accretion flow 
(ADAF; Narayan & Yi 1994, 1995; Abramowicz et al. 
1995) below certain critical accretion rate Mcr,ADAF and the 
luminous hot accretion flow (LEAF; Yuan 2001) above it 
(see review by Yuan & Narayan 2014). LHAF is thermally 
unstable (Yuan 2003). But if accretion rate is relatively 
low, Min ^ LHAF, the growth time-scale of the thermal 
instability is larger than the accretion time-scale, so the 
gas can remain hot throughout the flow. We call it type-I 
LHAF. Above Mcr.LUAF, the growth time-scale of the thermal 
instability is shorter than the accretion time-scale thus 
some cold dense clumps should be formed, embedded in the 
hot phase medium. We call it type-II LHAF or two-phase 
accretion flow. 

We follow the standard approach to calculate the 
dynamical structure and the emitted spectrum of ADAFs 
and Type-I LHAF (e.g., YCN05). The Lx at 3-9 keV band is 
then derived from the spectrum. For the two-phase accretion 
flow, optical/UV radiation from cold clumps will provide 
additional (likely dominant) seed photons for the inverse 
Compton scattering process, to generate the X-ray emission. 
The whole process is obviously complicated, depending on 
the detailed dynamics of the two-phases accretion flow such 
as the filling factor and temperature of clumps and so 
on. Following Yuan & Zdziarski (2004), we simplify the 
problem by replacing the electron energy equation with 
the Compton y-parameter. With the assumption that the 
electron advection term is zero in the two-phase flow, we 
can calculate the radiative cooling rate. The bolometric 
luminosity can then be derived (Xie & Yuan 2012). In order 
to derive 3-9 keV Lx Lx, we assume the spectrum to be 
a simple power-law with exponential cutoffs at both the 
high and low ends, i.e. Lg ~ exp(-Lmin/L) exp(-L/£n,ax), 
where the photon index T is constrained by the Compton 
y-parameter. is set arbitrarily to be 0.02 keV. £n,ax is 
numerically determined by the electron temperature, i.e. 
Lmax = LLe, at the location where most of the radiation comes 
out (Yuan & Zdziarski 2004). 


2.3 Basic parameters of the accretion—jet model 

In our calculations, the BH mass is fixed to Mbh = 10 Mq. 
In addition to the accretion rate Min of the accretion flow, 
our accretion-jet model contains several basic parameters as 
described below. 

Hot accretion flows are subject to strong wind (Yuan, 
Bu & Wu 2012; Yuan et al. 2015; Yuan & Narayan 2014), 
which makes the mass accretion rate in the accretion 
flow decreases with decreasing radius, following Min(r) = 
Min(f^out) (J^IRoutY, where index s is the wind parameter. 
Following numerical simulations of non-radiative accretion 
flow (see Yuan, Wu & Bu 2012 for the review of these 
simulations) we adopt s ~ 0.4 for the typical ADAF. With 
the increase of accretion rate, radiative cooling becomes 
more and more important. Consequently the wind strength 
may become weaker. We thus gradually reduce parameter y 
(down to 0.1) for high-M ADAF and Type-I LHAF. For the 
two-phase accretion flow we fix j = 0.1. 

The second is the viscosity parameter a. It is now widely 
accepted that magnetorotational instability is the mech¬ 
anism to transfer angular momentum (Balbus & Hawley 
1991, 1998). The value of a, however, is very diverse and 
not constrained. For example, Hawley et al. (2011) obtained 
a ~ 0.01 - 0.003, Penna et al. (2013) obtained a ~ 0.05 - 0.2, 
while Bai & Stone (2013) found that a can even be larger 
than unity. Likely the value of a mainly depends on the 
magnitude of the net flux of initial magnetic field (Pessah 
et al. 2007), i.e. a will be significantly larger when the net 
flux is nonzero. According to this result, the a value can be 
different in different sources, if the initial net magnetic flux 
contained in the accreting material from the companion star 
is different. 

The third is the plasma jj parameter, defined as the ratio 
of total pressure to the magnetic pressure. It constrains the 
strength of magnetic fields. We simply set ,6 = 10 throughout 
this paper, following numerical simulations of accretion flows 
(Yuan & Narayan 2014). 

The forth is the turbulent dissipation parameter 6, 
which characterizes the fraction of turbulent viscous heating 
that goes into electrons directly. Works have been done to es¬ 
timate the values of 6, by considering magnetic reconnection 
(Bisnovatyi-Kogan & Lovelace 1997; Quataert & Cruzinov 
1999; Ding et al. 2010), MHD turbulence (Quataert 1998; 
Blackman 1999; Lehe et al. 2009), or dissipation of pressure 
anisotropy (Sharma et al. 2007; Sironi & Narayan 2015). 
We now have the consensus that d ~ 0.1 - 0.5 but its exact 
value remains uncertain. From the observational side, by 
modelling the extremely dim supermassive BH, Sgr A*, 
Yuan et al. (2003) found that 5 ~ 0.5; while modelling to 
more luminous sources obtained a smaller value, d ~ 0.1 (Yu 
et al. 2011; Liu & Wu 2013). Since the sources we will deal 
with in this work are much more luminous than Sgr A*, we 
adopt d = 0.1. We will argue later in §3.1 that 5 = 0.1 is also 
favored in order to explain the radio/X-ray correlation. 

The fifth is a parameter that describes the coupling 
between the jet and the accretion flow. It is defined as the 
fraction of the mass accretion rate that goes into the jet, i.e., 


In literature, sometimes T/jd is assumed to be a constant. 
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Figure 1. The radio/X-ray correlation of BH X-ray binaries. As labelled in the figure, the points are observational data from three 
representatives. The three dotted lines show the fitting results to the data. The lines with different colour show the results calculated 
from the accretion-jet model. The (blue, green, and red) solid lines represent the pure hot solutions (ADAF and type-I LHAF), while the 
(green, magenta, and red) dashed lines represent the result of two-phase accretion model. The blue line corresponds to o- = 0.6, and others 
to O' = 0.06. The corresponding -Mn relationship is shown in Fig. 3. Note that the curves are superimposed at Lx^2 X lO^'^ergs^’. 


independent of Mi„(5f?s), perhaps for the reason of simplicity. 
Most MHD numerical simulations of jet formation have 
neglected radiation in hot accretion flow thus are scale- 
free to accretion rate. Physically, when the accretion rate 
is very low, radiation is not important and will not affect 
the dynamics of jet formation. In this case, a constant Tjjet 
is perhaps a good assumption. However, at higher accretion 
rates, radiation plays a more and more important role in the 
dynamics so rjj^t may no long be independent of the accretion 
rate. Given the uncertainties, we in our model set 77jet as a 
free parameter. By fitting the observed correlation, we find 
that 77jet decreases with increasing accretion rate (cf. Fig. 3. 
See also YCN05). This is consistent with the following two 
observational facts. One is that jet is suppressed in the soft 
state which has higher accretion rate compared with the 
hard state. Another is that the degree of radio-loudness of 
AGNs decreases systematically with increasing Eddington 
ratio (Fig. 10 in Ho 2008). 


3 MODELLING RESULTS 

3.1 General description to the method 

We note that, when Lx'^3x lO^^"* erg s“‘, both the ‘original’ 
and the ‘hybrid’ correlations have p ~ 0.6 (Coriat et al. 2011. 
Cf. Fig. 1). We first try to explain this part of correlation. 
The methodology is that we calculate the emission from 
the hot accretion flow, which dominates Lx, at various 
accretion rate Min and that from the jet, which dominates 
Lr, at different Mjet. The value of is adjusted to satisfy 
the correlation of this part. We emphasize that, the ~ 


M(5Rs) relationship is different for the two type of sources, 
whose viscosity parameters are different. The Pja ~ M(5f?s) 
relationship can be roughly described by power-law forms, as 
shown by the blue (for the ‘original’ correlation sources) and 
green (for the ‘hybrid’ correlation sources) solid lines in Fig. 
3. We then extrapolate these two power-law fitting functions 
to higher accretion rates and calculate the corresponding 
Lr and Lx to see whether we can explain the radio/X-ray 
correlations above Lx ~ 3 x 10’’’* erg s“* shown in Fig. 1. 

Before we introduce our results in detail, we first do 
some simple estimations, which is useful to understand the 
results presented in §3.2 and §3.3. We refer the readers to 
Fig. 1 of Xie & Yuan (2012). This figure shows the radiative 
efficiency of hot accretion flow as a function of the accretion 
rate. Because this figure is crucial to our present work, we 
reproduce this figure in Fig. 2, but with some modifications. 
One change is the value of a. In Xie & Yuan (2012), n = 0.1, 
now we have a = 0.06 and 0.6. The second is that now the 
wind parameter i is not a constant but is assumed to vary 
with the accretion rate. In addition, here we only adopt 6 - 
0 . 1 . 

Let us focus on the two curves in Xie & Yuan (2012) 
corresponding to h = 0.5 and 6 = 0.1. Each curve can 
be divided into three branches, with their boundary at 
M,,{5R,) ~ O.OlMEdd and M,,i5R,) ~ 0.004MEdd (Medd = 
10 LEdd/c'^ is the Eddington accretion rate)^. Both curves 


^ Note that Xie & Yuan (2012) used the net accretion rate M„et 
at the event horizon. Min(5Rs) ~ 1-9 Mnet, for the wind parameter 
s = 0.4 adopted by Xie &; Yuan (2012). 
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Figure 2. The radiative efBciency of hot accretion flow defined as 
E = Ctnt/Mi„(5Rs)c^, with Ltot being the total emission from the hot 
accretion flow. The red and blue curves are for a = 0.06 and 0.6, 
respectively. The red dashed curve is for the two-phase accretion 
flow. The radiative efficiency of SSD is also shown in this plot as 
grey dotted curve. 

can be described by a flat and very steep power-law forms 
at the high-Min and middle-Mjn branches, respectively. For 
the low-Mi„ branch, the d = 0.1 curve can also be roughly 
fitted by a power-law form, while the d = 0.5 curve is more 
curved. 

To explain the ‘original’ correlation like GX 339-4, 
which is described by a single power-law, it is natural to 
expect that the radiative efficiency of the hot accretion flow 
is also a single power-law, and the produced luminosity can 
cover the whole range of Lx, up to 4 x 10^’ ergs^' (refer to 
Fig. 1 of this paper). In this sense, the d = 0.1 branch is more 
promising than the 6 = 0.5 one, since the low-Mjn branch is a 
single power-law while the S = 0.5 branch is not. Of course, 
we also need that the low-A/i,, branch can produce Lx as high 
as 4 X 10’’’ ergs^'. 

To explain the ‘hybrid’ correlation like H1743-322, it is 
interesting to note that the three branches of the radiative 
efficiency curve may correspond to the three branches of 
the ‘hybrid’ correlations, namely the p = 0.6 correlation at 
low Lx, the transition regime at the middle Lx, and the 
p = 1.4 (or p = 0.98) correlation at high Lx, respectively. 
For example, the steep power-law efficiency curve means 
that with small increase of Mm, Lx will increase rapidly. 
This corresponds to the transition regime of the hybrid 
correlation. Again, the h = 0.1 curve in Xie & Yuan (2012) is 
better than the 6 = 0.5 one, since the steep power-law branch 
covers a wider range of radiative efficiency which is more 
suitable to explain the broad (about one order of magnitude 
in Lx) transition regime of the hybrid correlation. Below we 
present our detailed results. 

3.2 The ‘original’ correlation 

We adopt a = 0.6 and 6 = 0.1. The solid blue fine in Fig. 
1 shows the numerical result of the correlation (Note that 
below 2.5 X lO’’* ergs^*, the curves with different colors are 
superimposed with each other.). The observed ‘original’ p = 
0.6 correlation is reproduced, except at Lx> 10”ergs“* (but 
see our discussion at §4). The corresponding solid blue line 
in Figs. 2 and 3 show, respectively, the radiative efficiency 



Log[M,„(5R.)/M,,J 


Figure 3. The as a function of accretion rate MmiSRs). 
The colour and type of each line are the same to those in 
Fig. 1. Note that the red and green lines are superimposed at 
Min(5i?s)< 10-’“ MEdd. 

e and parameter ;7jet as a function of Mm(5Rs)- The typical 
value of !7jet(~ 10%) is consistent with what we have obtained 
in MHD numerical simulations of ‘disk-jet’ (e.g., Hawley & 
Krolik 2006; Ohsuga & Mineshige 2011; Yuan et al. 2015). 

The values of a is large compared with the value usually 
adopted in literature, but it is well within the reasonable 
range obtained from MHD numerical simulations (refer to 
§2). The reason for this choice is that, the critical accretion 
rate of ADAF, Mci.adaf ~ (0.1 - 0.3)a’MEdd (Yuan& 

Narayan 2014), will be large. Here the critical accretion 
rate measures the net accretion rate Min’- Therefore the 
luminosity from the ADAF can be high enough to roughly 
cover the full range of Lx of the original correlation. This, 
together with the fact that the radiative efficiency of an 
ADAF can be described by a single power-law in the case of 
d = 0.1 (cf. Fig. 2), explains why we can produce (nearly) a 
single power-law radio/X-ray correlation shown by the solid 
blue curve in Fig. 1. 

3.3 The ‘hybrid’ correlation 

We here choose a - 0.06 and d = 0.1 for this type of sources. 
We have tried several forms of zyjet - Mm(5Rs) functions, as 
shown by the (solid and dashed) red, green, and magenta 
lines in Fig. 3. The radio/X-ray correlation results are shown 
by the corresponding fines in Fig. 1. Below we will discuss 
these fitting results one by one. 

First, let’s look at the solid green fines in Figs. 1& 
3. The lines below and above ~ 2 x lO’^ ergs-' in Fig. 1 
correspond to the ADAF and type-I LHAF, respectively. In 
Fig. 3 the solid green line is similar to the solid blue fine 
(denoting the ‘original’ correlation) in the sense that both 
are straight lines (i.e., power-law function), but the radio/X- 
ray correlation results shown in Fig. 1 (solid blue and green 
lines) are different. The blue line is a single power-law while 
the green line becomes almost flat at Lx>2x lO’^ergs-'. 
This is because the dependence of radiative efficiency on 

’ The largest MmiSR^) in the solid blue line in Fig. 3 is slightly 
larger than 0.3«’MEdd- This is because wind parameter i adopted 
here is smaller than that in Yuan & Narayan (2014). 
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the accretion rate is different for ADAF and type-I LHAF, 
as we have explained in §3.1 (cf. Fig. 2). 

We have extrapolated the solid green line to higher 
Mi„(5Rs), i.e. the two-phase accretion mode, as shown by the 
dashed magenta line in Fig. 3. We then obtain a radio/X- 
ray correlation with correlation index p =» 1, as shown by the 
dashed magenta line in Fig. 1. This is in agreement with the 
correlation found by Gallo et al. (2012). 

However, theoretical considerations suggest that the Tyjet 
function may not extend as a simple power-law form as 
shown by the dashed magenta line in Fig. 3. This is because 
the accretion flow in the regime of dashed magenta line is 
in two-phase regime. For such kind of accretion flow, it is 
expected that, as the total (cold and hot) accretion rate 
increases, more fraction of the whole accreting gas will be in 
the cold phase. However, in this paper the accretion rate at 
the two-phase regime only takes the hot gas into account. 
Numerical MHD simulations of jet formation suggest that 
one important factor to determine the strength of jet is the 
magnitude of magnetic flux accumulated in the BH horizon 
or the innermost region of the accretion flow (§3.3 in Yuan & 
Narayan 2014), which is carried in by both the hot gas and 
the cold clumps. So for a two-phase accretion flow, with the 
increase of M„^(5Rs), Uja is more likely to be larger and larger, 
compared with the corresponding value shown by the dashed 
magenta line in Fig. 3. For simplicity we assume a constant 
Pjct, as shown by the dashed green line in Fig. 3. This results 
in a radio/X-ray correlation with p x 1.4, as shown by the 
green dashed curve in Fig. 1. This is consistent with that 
obtained by Coriat et al. (2011). 

The solid and dashed green curves can roughly repro¬ 
duce three branches of the ‘hybrid’ correlation with correct 
power-law index p, namely the low branch oi p x 0.6, the 
transition regime with p ~ 0, and the steep branch of p ~ 1.4. 
However, the fitting to the flat correlation branch, with 
Lx^ ~ 3 X 10^“* ergs^*, is clearly not satisfactory since the 
model overpredicts the Lr by a factor of ^2.5. 

To improve the fitting, we test another pjet- Mi„(5Rs) 

function, as shown by the (solid and dashed) red lines in 
Fig. 3. Compared with the green lines, the solid red line 
decreases faster with increasing Mi„{5Rs)- We also assume 
Pjet to be a constant at the two-phase accretion flow regime, 
i.e., Min(51?s)>0.02 The corresponding fitting results 

are shown by the (solid and dashed) red lines in Fig. 1. We 
see from Fig. 1 that now the model can explain both the 
p ~ 0 and the p = 1.4 correlations quantitatively well. 

The physical reason for the dashed red line in Fig. 3 
is the same with the dashed green line. Now the question is 
whether the solid red line physical? Or, given the uncertainty 
in jet formation theory, what information can we learn 
from this result? We note that the deviation of the solid 
red line from the solid green line in Fig. 3 begins from a 
relatively large Min close to Mct.adaf- This suggests that, as 
the accretion flow enters into the LHAF regime from the 
ADAF one, the radiative cooling becomes more important 
in affecting the dynamics of the accretion flow, including 
the jet formation process. This hypothesis can be tested in 
future by radiative MHD simulations of jet formation. 


4 SUMMARY AND DISCUSSION 

In term of the radio/X-ray correlation, Lr oc two 
types of sources have been observed. In the ‘original’ type, 
the correlation follows a simple correlation which can be 
described by a single power-law function, with the power- 
law index p fs 0.6. For the second type of sources, their 
correlation is ‘hybrid’, which is described by three power-law 
functions with different index p at the different luminosity 
regimes. In this paper, we try to explain these two types 
of correlation based on the coupled accretion-jet model of 
YCN05. This model has been successfully applied to model 
the multi-waveband spectrum of the hard state of BHBs 
and low-luminosity AGNs, and has explained the ‘original’ 
correlation (Yuan & Cui 2005; Yuan & Narayan 2014). 

In this model, the jet and the hot accretion flow are 
responsible for Lr and Lx, respectively. While Lr always 
follows a power-law function of the mass loss rate in the 
jet, the dependence of Lx on the accretion rate is much 
more complicated. Depending on the specific mode of the 
hot accretion, namely ADAF, type-I LHAF, and two-phase 
accretion flow, the relationship between the radiative effi¬ 
ciency and the accretion rate varies. This results in different 

values of p in the Lr -Lx correlation. The key assumption 

in this work is that, sources that always follow the ‘original’ 
correlation have a large viscosity parameter or; therefore 
the critical rate of ADAF is large and the accretion flow 
can always stay in the ADAF regime throughout the whole 
observed range of Lx. This explains why their correlation can 
be described by a single power-law form. On the other hand, 
the sources following the ‘hybrid’ correlation have a small a 
so the critical rate of ADAF is small. In this case, with the 
increase of accretion rate, their accretion modes change from 
ADAF to type-I LHAF, then to two-phase accretion flow. 
This explains why there are three branches of correlation 
for these sources. Our model indicates that the flat p ~ 0 
branch of the ‘hybrid’ correlation is actually X-ray bright 
(due to a quick enhancement in the radiative efficiency 
as Min increases) instead of radio faint. As suggested by 
MHD numerical simulations, the physical reason for the 
different value of a may be because of the difference in 
the net magnetic flux carried by the accreted material in 
the two types of sources. Unfortunately, the net magnetic 
flux of the accreting gas is unclear either observationally 
or theoretically. One possibility is that it may relate to the 
magnetic field configuration of the companion star. 

During the modelling, we find that the ratio between the 
mass-loss rate and the accretion rate pjet is not a constant 
of the accretion rate M,j,. Rather, in order to explain the 
observed radio/X-ray correlations, r/ja must decrease with 
increasing Mi„. Observationally, this is consistent with the 
facts that jets are present only in the hard state, and that the 
radio loudness of AGNs decreases with increasing Eddington 
rate. Theoretically, this indicates that radiation plays a 
certain role In affecting the jet formation. Detailed radiative 
MHD simulations of jet formation are required to test this 
hypothesis. 

Although a large a = 0.6 is adopted when we explain the 
‘original’ correlation, we still fail to reproduce the highest Lx 
(refer fig. 1). We note that GX 339-4 has almost the highest 
hard state Lx among all BH binaries detected. One comment 
is that the BH spin of this source may be large, a > 0.9 (e.g. 
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Miller et al. 2004; Reis et al. 2008; Yamada et al. 2009), 
while our calculations are for accretion on to a Schwarzschild 
BH. If a large BH spin were taken into account, we would 
produce a higher Lx since the radiative efficiency of accretion 
flow around spinning BH will be higher. 

One special source is Cyg X-1, where the BH accretes 
material from the wind of the high-mass companion. Its 
radio emission suffers additional free-free absorption by 
the stellar wind. After corrections of this absorption, the 
intrinsic radio/X-ray correlation of this source is likely 
p s! 1.4 (fig. 11 in Zdziarski et al. 2011). One caveat here 
is that some uncertainties exists in the correction of free- 
free absorption. If the p » 1.4 result is correct, under the 
scenario proposed in this work, it is likely that the viscosity 
parameter a of Cyg X-1 is moderately small, and the hard 
state of this source is described by the two-phase accretion. 
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